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This paper gives a pedagogical introduction to our knowledge of the spin structure 
of the nucleon. In particular, polarised deep inelastic lepton scattering is presented 
as a tool to study how the nucleon's constituents combine to generate its spin. The 
importance of semi-inclusive measurements is discussed and a window on future 
experiments in the field is given. 



1 Introduction 

This paper gives an overview of the spin structure of the nucleon with emphasis 
on a pedagogical presentation of the tools used to study this fundamental 
property of protons and neutrons. While many of the data currently available 
are shown, this paper is not intended to be an exhaustive summary of the field. 
Rather, the data are used to illustrate the concepts presented. The most recent 
progress in the field can be found in the. proceedings of the annual conference 
on deep inelastic scattering and QCDM. 

The last ten years have seen a wealth of measurements of the polarisation 
asymmetry in the cross-section for deep inelastic lepton scattering (DIS). This 
has been made possible by impressive advances in polarised target and beam 
technologies. Precise data now exist on inclusive polarised DIS from experi- 
ments at SLAC, CERN, and DESY. These data can be interpreted as showing 
that surprisingly, only about 30% of the nucleon's spin comes from the spins 
of the quarks. To make further progress on this problem, semi-inclusive mea- 
surements are needed where hadrons in coincidence with the scattered lepton 
are detected. These data allow us to determine the flavour of the quark that 
was struck and isolate the contribution to the nucleon's spin by valence and 
sea quarks. Furthermore, semi-inclusive data on charm production should give 
information on the contribution of gluons to the nucleon's spin. Current and 
future experiments are concentrating on semi-inclusive measurements. 

After a short discussion of magnetic moments, the formalism of polarised 



deep inelastic scattering is presented in Sec. 3.1-3.2. The data on the polarised 
structure function gi(x) are then shown in Sec. 3.3 which is followed by the 
interpretation of these data in terms of sum rules in Sec. ^. Section || shows how 
spin dependent structure functions can be used to determine the contribution 
of the quark spins to the spin of the nucleon. The Q 2 evolution of gi{x) 
is discussed in Sec. ||, while other spin structure functions are presented in 
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Sec. 0. A description of the experiments which have addressed the nucleon spin 
problem over the last decade is given in Sec. |8|. Semi-inclusive measurements, 
which are the current focus of the field are discussed in section |[ Finally, the 
future directions of the field are discussed in section [Ic]. 

2 The Spin of Three-Body States 

As a preamble to the discussion on polarised deep inelastic scattering, two 
simple examples of how three spin 1/2 particles combine to form a spin 1/2 
system are presented: the ground-state of 3 He and 3 H, and baryon magnetic 
moments. 



2.1 The Nuclear 3-Body Ground State 

The wave functions of 3 He and 3 H are among the most studied topics in physics. 
3 He is made up of two protons and one neutron, each of which has a spin of 
1/2. 3 H has one proton and two neutrons. The total wavefunction also has 
a spin of 1 /2 leading to the following possibilities for the combination of the 
three nucleons: 

L=0, S= 1/2 S-State (>90%) 

L=l, S= 1/2, 3/2 => P-State (<1%) 
L=2, S= 3/2 => D-State (8-9%). 

L is the orbital angular momentum of the 3-body system, while S is the total 
spin. Since the ground-state is overwhelmingly L=0, the two protons are anti- 
aligned due to the Pauli exclusion principle, and any spin effects observed 
in 3 He should come from the unpaired neutron. For example, consider the 
magnetic moment given by: 

V s =g s -s ; 5s(p)= 5.5855, 
g s (n)= -3.8327. 

We would have // s ( 3 Hc) = /i s (n) and /i s ( 3 H) = /it s (p) if the wave functions were 
purely S-state. This is approximately true in reality: 

Ms( 3 He) = -2.2174 w ^0)= -1.9130 
/i s ( 3 H) = 2.9788 w Ms(p)= 2.7927. 

A large part of the difference can be ascribed to meson exchange currents. 

Since these have opposite sign in the two nuclei, they should cancel in the 

average of the magnetic moments: 

|[M s ( 3 H)+/x s ( 3 He)]= 0.4257 
±[ Ms (p)+Ms(n)]= 0.4399. 
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Figure 1: Baryon magnetic moments. The proton, neutron, and A are used as input param- 
eters, while the other magnetic moments are predictions. The agreement of the data with 
the simple quark model (SQM) is very good. 

The remaining difference («3.2%) is due to the other components of the 
ground-state wavefunction (P- and D-states). 

As an experimental aside, the fact that polarised 3 He looks very much 
like a polarised neutron (at least as far as spin effects are concerned) makes 
it an ideal polarised neutron target. This has been exploited both at SLAC 
(E142/E154) and at DESY (HERMES). 

2.2 Baryon Magnetic Moments 

Since the simple model of three spin 1/2 particles combining to form a spin 
1/2 state works so well in the nuclear case, it is natural to apply it to the 
nucleoli. However, in this case we do not know the magnetic moments of the 
constituents, the quarks, since they cannot be isolated. Nevertheless, we can 
test the simple model by looking at the ratio of /i s (p) and /^ s (n) where the 
quark magnetic moments cancel. Experimentally, /x s (p)//i s (n)= -1.46, while 
the simple quark parton model predicts -1.50. Furthermore, we can use the 
proton, the neutron, and the A to set the values of fi s (u), /U s (d), /i s (s), and 
predict the magnetic moments of the other baryons. The results are shown 
in Fig. Ill The agreement between the predictions and the measurements is 
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Figure 2: Diagram of the Deep Inelastic Scattering process. The incoming lepton emits a 
virtual photon which is absorbed by one of the quarks in the nucleon. 



very good indicating that the simple model works well. However, while the 
magnetic moments are consistent with the model presented here, they depend 
only on the static properties of the nucleon wavefunction (more precisely at 
Q 2 = 0) and the rest of this paper will discuss how this picture is too simple. 

3 Inclusive Polarised Deep Inelastic Scattering 

3.1 Polarised DIS Formalism 

The main process we will be concerned with in these lectures, deep inelastic 
scattering, is depicted in Fig. |^. An incoming lepton (shown here as a positron 
or an electron) emits a virtual photon which is absorbed by a quark in the 
nucleon. The nucleon is broken up and the struck quark and the target remnant 
fragment into hadrons in the final state. Only the lepton is detected in inclusive 
measurements while detection of one or more hadrons in the final state (semi- 
inclusive data) adds important information on the scattering process. W ± and 
Z° exchanges are not important at the energies of the current spin experiments. 
The kinematic variables relevant for this process are listed in Table [lj The 
formalism for DxjcijL Loplastic Scattering (DIS) is developped in many texts on 
particle physics fflMfl. 

The DIS cross-section can be written as follows: 

d 2 a a 2 E' 
= - A -L, V W^ (1) 



dfl dE' 2M p Q 4 E 



4 



Table 1: Kinematic variables relevant to Deep Inelastic Scattering 



(k,E) 


Momentum and energy of the incoming lepton 


(k',E>) 


Momentum and energy of the outgoing lepton 


v = E-E' 


Energy transfer in the reaction 


q = k-k' 


Four momentum transfer 


Q 2 = -Q 2 


= AEE' sin 2 (0/2); 6: e ± scattering angle 


x = Q 2 /2M p v 


Fraction of the nucleon momentum carried by 


the struck quark 


y = v/E 


Fraction of the incident energy transferred 


W 2 = M 2 V + 2M p v - Q 2 


Square of the invariant mass of the 


photon-nuclcon system 


z = E h /v 


Fraction of the energy transferred carried by 




a particular hadron h 


x F = 2p n /W 


Scaled longitudinal momentum of a hadron 



where a is the fine structure constant and M p is the mass of the proton. No- 
tice that the cross-section falls off very quickly as a function of momentum 
transfer Q. is a tensor which describes the emission of the virtual photon 
by the lepton; it can be calculated exactly in QED. W^ u is a tensor describing 
the absorption of the virtual photon by the target; it contains all of the infor- 
mation related to the structure of the target. Symmetry considerations and 
conservation laws determine the form of W^ v : 

W» v = -g» v Fi + j/y/i/.Fa (2) 
+ %/v e^ XrT q X s a ■ gl + ijv 2 e^ Xa q X (p.q s a - s.q p a ) ■ g 2 . 

Fi and F 2 are unpolarised structure functions, while gi and g 2 are polarised 
structure functions which contribute to the cross-section only if both the target 
and the beam are polarised. In general, the structure functions depend on v 
and Q 2 . However, if the scattering occurs from pointlikc constituents, they 
depend only on one variable x = Q 2 /2M p v. This is the phenomenon of scaling 
which was the first experimental confirmation that the nucleon is composed 
of smaller particles. In the infinite momentum frame, x is shown to be the 
fraction of the nucleon's momentum carried by the struck quark. The structure 
functions can be interpreted in the quark-parton model in a straightforward 
way: 

= \ E4 (tfo*) = \ E4 «/(*) ( 3 ) 

/ / 
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Figure 3: Diagram of Polarised Deep Inelastic Scattering. The virtual photon is always 
absorbed by a quark with opposite spin direction. The top panel shows the process that is 
sensitive to quarks polarised parallel to the nucleon's spin (q~f(x)), while the bottom panel 
is for quarks with their spin direction anti-parallel to the nucleon (qj(x)). 

9i{*) = \ E e / = \ E4 A ?/W ( 4 ) 

/ / 

where the sum is over quark flavours and e/ is the charge of quark, qf (qj ) 
is the distribution of quarks and anti-quarks with their spins in the same 
(opposite) direction as the spin of the nucleon. F± (x) measures the momentum 
distribution of the quarks in the nucleon, while g\ (x) gives the distribution of 
the spins of the quarks. For spin 1/2 quarks, F 2 is related to i*\ by El 

F 2 {x) = 2xF 1 {x). (5) 

Finally, g 2 — in the quark-parton model since it is related to transverse 
degrees of freedom which are absent in the simple QPM. 

An intuitive picture of how gi(x) is sensitive to spin degrees of freedom 
and how to measure it is given in Fig. [|. The incoming polarised lepton at the 
left emits a circularly polarised virtual photon with spin projection pointing 
to the right in this picture. This photon can be absorbed only by a quark 
with spin projection in the opposite direction because the final state (a quark) 
must have spin 1/2 and can therefore not have spin projection 3/2. In the top 
panel of Fig. the target polarisation 'points' to the left so that we measure 
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Figure 4: Diagram of the angles relevant to Polarised Deep Inelastic Scattering. See the text 

for details. 



the cross-section for beam and target polarisations anti-parallel. This cross- 
section is labelled cr^ and is sensitive to the distribution of quarks with their 
spins in the same direction as the spin of the host nucleon, or q~j(x). In the 
bottom panel of Fig. |], the target spin direction has been reversed so that 
we measure the cross-section for beam and target polarisations parallel (ct^). 
However, the elementary process is not changed. The photon absorbs on a 
quark with spin projection 'pointing' to the left and this process is sensitive to 
qj{x). Since gi(x) is related to the difference of q~f(x) and qj(x) (see Eqn. 
it can be deduced by measuring the polarisation asymmetry: 



J1 



-4„ 



(6) 



Further insight can be gained into how to measure gi(x) and 52 (x) by 
considering Fig. [|. k and k' are the usual lepton momenta, <fi is the angle 
between the lepton scattering plane and the target spin, and a is the angle 
between the beam and target spins. The target polarisation is said to be 
longitudinal (transverse) when a = 0° or 180° (= 90° or 270°). The difference 
in the cross-sections (Act) for the two relative directions of the beam and target 
spins can now be written as follows: 



d 3 Aa 
dxdyd(f> 



cosa • {a gi(x) + b 52(2;)} 
cos <j> sin a ■ c ■ { | g\ (x) + gi (x) } . 



(7) 



Each of the two terms in braces can be isolated by preparing the target with 
either longitudinal or transverse target polarisation. The kinematic constants 
a, b, and c are of order 0(1), 0(1/Q 2 ), and 0(1/Q) respectively. Therefore, 
g\(x) dominates for longitudinal target polarisation, even at moderate Q 2 . 
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This statement is strengthened by the fact that there are good arguments that 
g2{x) is small since it is related to transverse degrees of freedom in the miclcon. 
The small value of 52(2;) has been confirmed by experiment as we will see in 
Sec. Once gi(x) has been determined, transverse target polarisation can 
be used to maximize the sensitivity to #2 (%)■ However, notice that in this case 
the asymmetry is reduced by the factor c oc 1/Q, making this measurement 
more difficult. 



3.2 What is measured in practice? 

This section gives the steps required to go from a measurement of the cross- 
section asymmetry to the polarised structure functions. First, one must realize 
that the experimentally observed asymmetry (A meas ) is reduced because we 
do not have perfect beam and target polarisation. It is related to the true 
asymmetry (^4y) by: 

Ameas = jyU Jfft = (/d ' PB ' Pt) ' A\\ , (8) 

where pb and pr are respectively the beam and target polarisations, and fn is 
the target dilution factor. fr> is defined as the fraction of polarisable nucleons in 
the target (1 for H and D gas; \ for 3 He gas; f° r N H 3 ; and |2 for C4H9OH). 
fo is further reduced by extraneous materials in the target such as windows, or 
vapour introduced for the optical pumping process. Beam polarisations vary 
from 0.5 to 0.85, while the range of available target polarisations is 0.4-0.9. 
The ratio of the measured to the real asymmetry is in the range 0.07-0.12 
(CERN); 0.17-0.35 (SLAC); 0.25-0.50 (DESY) (see Sec. |). Since some of the 
asymmetries are small to begin with, the experiments must be very accurate, 
requiring high statistics. 

The measured asymmetries for longitudinal and transverse target polari- 
sation, corrected for the factors in Eqn. |[ are written in terms of the cross- 
sections as follows: 

a U — o-TT ortJ. — o-TT f 1 da 

A u = a n + (J ii ; Al - = an + £tt ■ where b = / ^ 



COS ( 



While Am and A± depend on beam and target related polarised cross-sections, 
the physically significant asymmetries for photon absorption on the nucleon 
level are given by: 



(T 1 —(7 3 



A 1 = ^—^ ; A 2 = - n = — > ( 10 ) 

(Ti + (73 ^[Ci + C73 ) Gt 



where o~i is the absorption cross-section when the photon and the nucleon 
have their spins anti-aligned, and (J3 the cross-section when they are aligned, 
or and <7tl are the absorption cross-section for tranversely polarised virtual 
photons and the longitudinal-transverse interference respectively. 

The measured asymmetries are related to the photon absorption asymme- 
tries by the following equations: 

A\\=D-(A 1 +nA 2 ) ■ A ± =d-{A 2 -$At). (11) 

D, r/, d, and £ are depolarisation factors given by: 



y 2 +2(1 -y)(l + R) 2-y M 1 + e 



S, = "h—ft— ; e = — - — ; i = Ql v = 2mx/Q. 

2e 1 - y + TL. 

R is the ratio of longitudinal to transverse cross-sections I. The depolarisation 
factors account for the loss of polarisation from the lepton to the virtual photon 
and depend on kinematics. 

Finally, the photon asymmetries are related to the structure functions as 
follows: 

A 91^92 A 7(ffi+g 2 ) (13) 

Ft ' Fi K ' 

If gi and A 2 are ignored, the equation for g\ simplifies considerably: 

91 ( X ) n Fl (x) ■ At(x) * Ft(x) ■ (14) 

3.3 Data on the Polarised Structure Function g\{x) 

The EMC results^ published in 1988, which indicated that little of the spin of 
the proton comes from the spins of the quarks, generated a series of experiments 
at CERN (SMC: H & D targets), SLAC (E142 & 154: 3 He, and E143 & 155: 
H & D), and DESY (HERMES: H, D, & 3 He) to study the "nucleon spin 
puzzle". As a result, we now have precise data on gi(x) over a range of Q 2 . A 
selection of these data, as of the time of this workshop (Feb. '98), is shown in 
Figs. (see the figure captions for references). There are several choices for 
the form of these plots. Here, the abscissa is plotted on a logarithmic scale to 
allow the points at low-x to be seen clearly. xg\ (x) rather than g\ (x) is plotted 
on the ordinate because this is the appropriate quantity to plot to illustrate 
the contribution of each data point to the integral of g\{x), the importance of 
which in terms of the spin content of the nucleon will be discussed in Sec. ||-|[ 
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Figure 5: Data on xg\ (x) from SMC 11 , E143 14 , and preliminary results from HERMES 16 . 
The E143 data arc at Q 2 = 3 GeV 2 , while the SMC and HERMES data arc at the measured 

Q 2 for each x bin. 
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Figure 6: Data on xg^(x) from SMC 12 and E143 15 . The E143 data are at Q 2 = 3 GeV 2 , 
while the SMC data are at the measured Q 2 for each x bin. 
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Figure 7: Data on xg^(x) from E143 15 , HERMES 17 , and E154 18 . The E143 data are at 
Q 2 = 3 GeV 2 , E154 is at Q 2 = 5 GeV 2 , and the HERMES data are at the measured Q 2 in 

each x bin. 



A few general points can be made about the data. First, although they 
are taken with widely varying experimental techniques (see Sec. the data 
are consistent with each other, giving confidence in the results. Second, while 
the Q 2 varies substantially for the different experiments ((<3 2 ) = 2-10 GeV 2 ), 
gi(x) shows relatively little Q 2 -dependence over most of the x-range, at least 
within the experimental uncertainties. Third, while the SMC data have worse 
statistical precision than either the SLAC experiments or HERMES, they ex- 
tend to lower x which is important for the extrapolation of gi{x) to x = 0. 
Finally, the precision of the data is quite good and the experiments are now 
close to being limited by systematic uncertainties. 

One should notice that while g\(x) is large and positive for 0.1 < x < 0.6, 
gi(x) is small and negative for all x. The reader is referred to the experimental 
papers for a more complete presentation of the data. 



4 The Interpretation of 51(2;) in Terms of Sum Rules 

The spin structure functions are often interpreted using sum rules which relate 
the integrals of g\(x) over all x to the polarisation of the quarks in the nucleon. 
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4-1 The Bjorken Sum Rule 

If we use Eqn. |] for the proton and the neutron, we get: 



2 V 9 



-Ad- 



r?= f gl{x)dx 
Jo 

= f dx=\ (~Au + ^Ad + ±As\ , 

where we have used isospin symmetry to make the simplifications: 

Au p = Ad n — Au 



Ad p 
As z 



Au n 

■ As r , 



Ad 
As 



(15) 
(16) 



(17) 
(18) 
(19) 



with the corresponding equations for the anti-quarks. Au, Ad, and As in 
Eqns. [l6| are therefore for the proton. The sea cancels in the difference of 
the two integrals so that 



[gV(x)-g?(x)]dx = ±(Au-Ad) = ± 



9A 



gv 



(20) 



where the final step comes from neutron /3-decay, and <?a and gv are the axial 
vector and vector coupling constants. This sum rule, first derived by Bjorkenli3 
using current algebra, has little model dependence and is fundamental to QCD. 
If higher order QCD effects are_taken into account, the full expression for the 
Bjorken Sum Rule is obtainedl23: 



[gl{x)-gUx)]dx = - 

D 



fjA 



Corr 



Corr = 1 — 



a s (Q 2 ) 



-20.2153 



a s (Q 



7T 

2-.\ 3 



- 3.5833 



0(130) 



gv 

7T 

« S (Q 2 ) 



(21) 



where a s is the Q 2 dependent strong coupling constant. The precise data on 
gi(x), gi(x), and gf(x) can be used to test the Bjorken sum rule. The current 
situation is shown in Fig. || which is a plot of the neutron integral (r™) versus 
the proton integral (1^). The overlap of the bands for the three targets is 
consistent with the sum rule to one standard deviation. The Bjorken sum rule 
has been verified to about 10%. 
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E143 

Figure 8: Plot of r™ versus T^. The predictions of the Bjorken and Ellis-Jaffe sum rules are 
shown on the diagonal band from the lower left to the upper right of the figure. While the 
data and the Bjorken sum rule overlap within one sigma, the Ellis-Jaffe prediction is roughly 
two sigma away from the overlap region in the data. 



4-2 The Ellis-Jaffe Sum Rules 

Because of the difficulty of producing a polarised neutron target, it became 
clear that data on the ucpton would be available much before data on the 
neutron. Ellis and JaffeEil derived sum rules for rf and T" separately using 
SU(3)f flavour symmetry and the assumption that As= 0. One can write the 
following from neutron and hyperon /3-decays: 

Au - Ad = F + D (22) 
Am + Ad - 2As = 3F - D 



where F and D are SU(3)p coupling constants. If As = 0, we can solve Eqs. 22 
for Am and Ad and write Eqs. [l5], |l6| in terms of F and D: 

f gl{x)dx = ^F-D) (23) 
J g?(x)dx=±(6F-4D) (24) 
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Table 2: Values of the Ellis-Jaffe integrals from the various polarised DIS experiments. The 
second number in the fourth column is the statistical uncertainty on the measurement, while 
the third number is the systematic error. 



Experiment 


Integral 


(Q 2 ) (GeV^) 


Measurement I 




f. 


EMC 




10.7 


0.114 ± 0.012 ± 0.026 


ii 




SMC 


r? 


10 


0.136 ± 0.013 ± 0.011 


ii 






rf 


10 


0.041 ± 0.006 ± 0.005 


12 






1 1 


10 


-0.046 ± 0.018 ± 0.019 


12 




E142 


1 1 


2 


-0.031 ± 0.006 ± 0.009 


13 




E143 


r? 


3 


0.127 ± 0.004 ± 0.009 


U 






rf 


3 


0.042 ± 0.003 ± 0.004 


IS 






pn 


3 


-0.037 ± 0.008 ± 0.011 


L5 




E154 


1 1 


5 


-0.041 ± 0.004 ± 0.006 


18 




HERMES 


rn 
1 1 


2.5 


-0.037 ± 0.013 ± 0.008 


17 





These sum rules are usually written in a different form: 



gl{x) dx= — 



fjA 



9v 
9a 



D 



9v 



5 3F 

H 

3 F + D 



5 3F-D 
1 + 3 F + D 



x QCD Corr. 
x QCD Corr. 



(25) 
(26) 



Using the fact that F + D = gA/gv — 1.257 ± 0.003 from neutron /3-decay 
and the determination of Fj D from a global fit to hyperon /3-decays (F j D = 
0.575 ± 0.0161 we can get numerical predictions for the Ellis-Jaffe integrals. 
For exampleEj, at Q 2 = 5 GeV 2 : 

rf (EJ) = 0.172 ± 0.009 ; T1(EJ) = -0.018 ± 0.009 

where higher order QCD corrections are included. While the Bjorken sum rule 
is verified by data, there is significant disagreement between the measured val- 
ues of r? and T" and the Ellis-Jaffe predictions. The situation is summarized 
in Fig. |9j where deviations of several a can be seen between the data and the 
predictions. The values of the Ellis-Jaffe integrals for the different experiments 
are given in Table 2. 

The disagreement of the data with the predictions of Ellis and Jaffe is not 
surprising for several reasons. From the theory standpoint, the assumptions of 
SU(3)f symmetry and As = are crude at best. On the experimental side, 
there is great controversy over how the extrapolation to the unmeasured region 
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Figure 9: Results for rj, T™, and Ff compared to the predictions of Ellis and Jaffe (shaded 

area) . 
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at low x should be done. Regge theory has been used in the past but recent 
low x data (F2 and g\) have shown this to be inappropriate. Several other 
functions have been tried in the low x region butihere is no agreement to date 
on the proper form to use (see for example Ref.E3 and references therein). 



5 The Spin Content of the Nucleon 

The expressions for the integrals T\ and T", and for the matrix elements of 
hyperon /3-decays in terms of Au, Ad, and As can be used to deduce the 
fraction of the nucleon spin due to the spins of the quarks. In practice, global 
fits are made to all the data. However as an illustration, a simple example is 
shown here which has the added advantage of giving_historical insight into the 
'Spin Crisis' since it uses the original EMC data El Consider the following 
system of three equations and three unknowns: 

T{ = -{-Au+-Ad+-As) = 0.114 ±0.029 (EMC) 
2 9 9 9 

Am - Ad = 1.257 ± 0.011 (neutron j3 - decay) (27) 
Am + Ad - 2 As = 0.579 ± 0.016 (hyperon (3 - decay). 

Solving for Au, Ad, and As, we get: 

Au = 0.705 ±0.09 ; Ad = -0.553 ± 0.09 ; As = -0.214 ± 0.09. (28) 

The total contribution to the nucleon's spin from the quark spins is therefore 

AS = Au + Ad + As = -0.062 ± 0.16, (29) 

consistent with and several standard deviations from what was expected. 
The EMC did a somewhat different analysis (assuming the validity of the 
Bjorken sum rule) but they also obtained AS w 0. This surprising result lead 
to the new round of experiments designed specifically for polarised DIS. Now 
that many more data exist, global fits are done to get the Ag's with resulting 
values of AE in the range 0.2-0.3, still significantly below expectation. 

In the naive quark parton model, all the spin of the nucleon comes from 
the spin of the quarks so that AE= 1. If relativistic wave functions are used 



x 

<7-p 



E+m ' X 

the lower components contribute orbital angular momentum and AE is reduced 
to 0.75B. The Ellis- Jaffe prediction for AE can be obtained from Eqn. |22| with 
As=0: 

Au + Ad = 3F-D =s> AE = 3 x 0.459 - 0.791 = 0.586. 
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It is clear that no matter what one does, the measurements of AS are all 
in disagreement with the predictions indicating that there is still something we 
do not understand about the spin of the nucleon; or that there is something 
wrong with our interpretation of AS (see below) . 

6 The Q 2 Evolution of the Structure Functions 

6.1 Alternate expression for g\{x) 

Eqns. |], [l^, and 16 give simple expressions for g\(x) in term of the quark 



distributions. These can be generalized to include gluons and written_in a 
form which is convenient for the Q 2 evolution of the structure functions EJ: 

(e 2 ) 

gx(x, Q 2 ) = [C ns ® Aq ns + C s <g> AS + 2n f C g ® Ag] , (30) 



11 f 

where (e 2 ) = — (31) 

' i=l 

Uf ( e 2 \ 

Aq ns = ^J^^-lJ (Aqi+Aqi) (32) 
3 

= -(Am — Ad) (for 2 flavours), and 
5 

AS = ^(Aft + Ag,)- (33) 
i=l 

Aq ns is called the non-singlet quark distribution, AS is the singlet quark dis- 
tribution, Ag is the gluon polarisation, nf is the number of quark flavours, and 
Cns, C s , and C g are coefficient functions calculable in QCD. The ® in Eqn. 30 
refers to a convolution defined as follows: 



{C®q){x,Q 2 )= f c[-,a s ) q(z,Q 2 )^. (34) 



The advantages of this formulation can be seen if one considers thj^ evolu- 
tion of g\ with Q 2 . This evolution is given by the DGLAP equations cB: 

±Aq ns = ^Lp™®Aq ns ; t= log (35) 
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Pqq and the P?j are splitting functions and A is the constant of the strong 
interaction (w 200-300 MeV). The first thing to notice is that the Q 2 evolution 
of the non-singlet term does not depend on the gluons. Since the Bjorken sum 
rule is related to the non-singlet term (see Eqns. ^C], [52]), it depends only on 
the quark distributions. On the other hand, the evolution of the singlet term 
depends on both the quark and gluon distributions. Therefore, we must be 
very careful in our interpretation of AS. In fact, the quark distributions as 
measured in DIS (Aq) are modified by a gluon term A(/c3: 

Aq = Aq-^Ag. (37) 

nf is the number of active quark flavours (3 if heavy quarks are ignored) . This 
is thought to be the best candidate to explain the relatively small measured 
value of AE, and the search for an indication that the gluons are polarised has 
become a focus of the field. 

6.2 Gluon Polarisation from Q 2 Evolution 

There have been attempts to extract a value of Ag from the Q 2 evolution of 
gi(x) by using Eqns. [SO] [35] and |6]to perform a global fit to data at different 
Q 2 . The result of these fits for the E143 dataEZl is shown in Fig. |l0] where some 
Q 2 dependence is seen at low Q 2 , albeit jSmall-J^alues of Ag in the range 1.5- 

2.0 have been deduced in several studiesl25E§Eal23. However, the uncertainties 
are large and the data are consistent with Ag=0. Much better data will be 
required for this technique to produce reliable results. 

7 Other Structure Functions 

7.1 g 2 (x) 

Until now the discussion has centered on g\. However, several other structure 
functions can be measured. g 2 {x) has already been introduced in Eqn. |] and 
Sec. |3.2[ While g 2 could be important as a correction to g\ , it is also interesting 
in its owu-right because it could offer a way of isolating the effects of a twist-3 
operator Ell. Twist is a number that arises in the operatpc product expansion 
and is an indication of the order of the operator involvedEa. Twist-2 operators 
are leading-order, while twist-3 are next-to-leading order, etc. The following 
expression for g 2 shows a term Q72) which is pure twist-3: 

92(x, Q 2 ) = g¥ W (x, Q 2 ) + h(x, Q 2 ) (38) 
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Figure 10: gf/F^ as a function of Q 2 in separate x bins. The curves are different global fits 
to the asymmetries. This figure is taken from Ref. 27 where more detail can be found on the 

fits. 
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Figure 11: Measurements of (a) xg% and (b) xg% from E143 at SLAC 34 . Systematic un- 
certainties are indicated by the bands. The solid curve shows the g^ w calculations. Bag 
model calculations at Q 2 = 5.0 GeV 2 by Stratmann 35 (dotted), and Song and McCarthy 36 

(dashed) are also shown. 



with gY W (x,Q 2 ) = ~ 9l ( X ,Q 2 )+ f 9 -±^ldy (39) 

Jx y 

The term g^ w was derived by Wandzura and WilczekEl. It is twist- 2 (i.e. 
leading order) and can be calculated from g\. A deviation of g 2 ^ om r-SiY W 
would signal the effect of a twist-3 operator. Data on g 2 from SLAC a are 
shown in Fig. |ll|. The data are consistent with either or g\ vw which is 
shown as the solid line. They are not precise enough to distinguish between 
these two possibilities or to isolate any possible twist-3 contribution. However, 
these data confirm that corrections to gi due to g 2 are small. More precise 
data on g 2 are expected from SLAC in 1999. 
Burkhart and CottinghamE3 have derived a sum rule for g 2 : 

/ g 2 (x)dx = Q (Q 2 ~>oo). (40) 
Jo 
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The data in Fig. |ll| have been integrated with the following results 0: 

/ g p 2 {x) dx = -0.013 ±0.028 (41) 

JO. 03 

f g$(x) dx = -0.033 ± 0.082 (42) 

Jo. 03 

both consistent with 0. 

7.2 Still more Structure Functions 
hi(x) 

hx(x) is a twist- 2 structure function like F\(x) and <7i(x)il. It is related to 
tensor currents in the nucleon. However, it is more difficult to measure because 
it is chiral-odd. A longitudinally polarised beam and a transversely polarised 
target are needed to measure h\{x). Furthermore, because this structure func- 
tion is chiral-odd, pions must be detected in coincidence with the lepton to 
conserve chirality. 

h(x) 

b\(x) is measured with an unpolarised beam and a longitudinaly polarised D 
target in each of three substates (md= 1, 0, -1). It is given by the expression: 

& i (*) = \ E 4 (*) + «7 (*) - H (*)] (43) 
/ 

and measures quarks from nuclear binding El. 
A(x) 

A (a;) is also measured with an unpolarised beam and a D target, but in this case 
the target is transversely polarised in the md—0 state. The cp distribution of the 
scattered leptons determines A(x) which is sensitive to gluons not identified 
with individual nucleonsc3. 



8 Experiments to Study the Spin Structure of the Nucleon 

We will now spend some time describing the recent experiments which have 
studied polarised deep inelastic scattering. The emphasis will be on the unique 
features of each one. 
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Figure 12: SMC Spectrometer. 



8.1 CERN: The Spin Muon Collaboration (SMC) 

In the SMC experiment, 100-190 GeV muons scatter from butanol (C^^Orl), 
deuterated butanol (C4D9OD), and ammonia (NH 3 ) targets (see Ref. liil and 
references therein). The experiment is an extension of the pioneering work 
done at CERN in the eighties and uses essentially the same spectrometer as 
the EMC, shown in Fig. [l2]. This device is centered around a large aperture 
dipole magnet which provides a 2.3 T.m magnetic field. Particle tracks are 
determined by a series of multi-wire proportional chambers and drift cham- 
bers. Muon identification is done using a thick hadron absorber behind the 
tracking chambers, followed by streamer tubes. The trigger is generated by 
three hodoscopes behind the absorber. 

The muon beam is produced by the decay in flight of pions and kaons. 
The current is equivalent to about 0.5 pA (4.5 x 10 7 muons per spill in 2.4 s 
every 14.4 s). The muons are polarised because of parity violation in the 
pion and kaon decays. The typical beam polarisation was 80%, determined by 
a calculation and confirmed by a measurement of the Michel spectrum of /1 
decays downstream of the experiment. 

The target is polarised using dynamic nuclear polarisation. In this tech- 
nique, the high polarisation of electrons in paramagnetic impurities is trans- 
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Table 3: Running conditions for SMC. 



Year 


E B (GeV) 


PB 


Pt 


Id 


Target 


'92 


100 


0.84 


0.35 


0.23 


d-butanol 


'93 


190 


0.79 


0.86 


0.13 


p-butanol 


'94 


190 


0.79 


0.50 


0.23 


d-butanol 


'95 


190 


0.79 


0.50 


0.23 


d-butanol 


'96 


190 


0.79 


0.89 


0.17 


p-NH 3 



ferred to protons using simultaneous flips of these spins induced by microwave 
irradiation. The polarisation direction is chosen by adjusting the frequency 
of the microwaves. The target operates in a high magnetic field (2.5 T) at 
a temperature of 50 mK. The target orientation can be chosen longitudinal 
or transverse to the incoming muon. One of the remarkable features of the 
experiment is that two targets with opposite spin direction are used simulta- 
neously. Each target cell is 60 cm long and 5 cm in diameter. The use of two 
targets is ideal for measuring yield ratios and asymmetries, greatly reducing 
the systematic uncertainties related to luminosity and detector performance. 
On the other hand, the SMC targets have the disadvantage of containing a 
large fraction of unpolarised nucleons which greatly dilutes the experimentally 
observed asymmetry. A measure of this effect is given by the dilution factor 
fr> which is the ratio of polarisable nucleons divided by the total number of 
nucleons in the target. The factor fn varies from 0.13 to 0.23 in the SMC 
experiment. The polarisation of each half of the target is flipped from time to 
time to take into account the difference in the acceptance of the spectrometer 
for each target section. The polarisation is measured using NMR and was 
found to be 86% for protons and 25-50% for deuterons. Target thicknesses 
of about 5 . 10 24 cm -2 were obtained leading to luminosities on the order of 
1-2 . 10 31 cm~ 2 s _1 . The running conditions for SMC from 1992 to 1996 are 
summarised in Table ||. 

Another unique feature of SMC is that it has the highest beam energy of 
the polarised DIS experiments and therefore the highest Q 2 ((Q 2 ) ~ 10 GeV 2 ) 
and lowest x {x m i n w 0.003). The low x points are especially valuable in 
determining the best functional form to use for the extrapolation of g\{x) to 
x = 0. 
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8.2 SLAC: E142 - E143 - E154 - E155 



The SLAC experiments use spectrometers in end-station A. While the technol- 
ogy used for the polarised electron beam is the .same for all these experiments, 
the targets are very different (see references inE3 for more details). 

The beam is produced by photoemission from an AlGaAs crystal induced 
by circularly polarised laser light. The direction of the electron spin is chosen 
at random on a pulse by pulse basis. The beam polarisation is measured by 
M0ller scattering from thin ferromagnetic foils. It was about 36% for E142 
and improved dramatically for the subsequent experiments to 82-86% by using 
a strained AlGaAs crystal. The beam energy varied from 9.7-29.1 GeV for 
E142-143 to 48.3 GeV for E154. Beam currents also varied from 0.5-3.5 fiA 
(e.g. 2.10 11 e" per pulse with a 120 Hz repetition rate for E142, to 3-9. 10 10 
e~ per pulse for E154). 

E142 and E154 used a 3 He target 0, while E143 and E155 used NH 3 , 



ND3 and LiD targets. As discussed in Sec. 2.1, any spin effects observed 
in 3 He scattering can be ascribed to the neutron. The SLAC target uses 
optical pumping by infrared lasers of Rb atoms, followed by spin-exchange 
with the 3 He nuclei. Two cells are used, one for optical pumping and one 
for the scattering, connected by a thin tube. This helps reduce the amount 
of Rb in the beam which in turn keeps the dilution factor fu relatively high 
{Id ~ 0.33 for E142). The target cell was 30 cm long by 2 cm in diameter and 
contained between 8.4 atm (E142) and 9.5 atm (E154) of gas at 20°C. The 
polarisation, measured using NMR, was 33% for E142 and 38% for E154. The 
target polarisation was reversed from time to time to reduce the systematic 
uncertainty in the asymmetry. The NH3 and ND3 targets used by E143 and 
E155 (LiD instead of ND3) are similar to the SMC target in that they also use 
dynamic nuclear polarisation. However, in this case only one target cell is used 
and it is much smaller than the SMC target (3 cm long and 2.5 cm in diameter). 
The target field is 4.8 T and the temperature is 1.6 K. The polarisations were 
in the range 65-80% for the proton and «25% for the deuteron. The dilution 
factors range from 0.17 (H) to 0.25 (D). The advantage for E155 of LiD over 
ND3 comes from the fact that the Li, which can be described approximately 
as an a-particle plus a deuteron, is also polarised. Hence the dilution factor is 
better for LiD. 

The spectrometers used for E142-143 are shown schematically in Fig.|l3|. 
The electrons are detected at scattering angles of 4.5° and 7.0°. The spectrom- 
eters were modified for operation at higher energy for E154-155. In particular, 
the scattering angles were changed to 2.75° and 5.5° and a new 10° spectrom- 
eter was built. Electrons are identified by Cerenkov detectors and lead-glass 
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Figure 13: Spectrometer for SLAC Experiments E142 and E143. 

calorimeters while tracking is done with scintillator hodoscopes. The geo- 
metrical acceptance of these spectrometers is small compared to the other 
experiments. 

The running conditions for the SLAC experiments are summarized in Ta- 
ble 1 



8.3 DESY: HERMES 

The HERMES experiment uses a polarised positron (electron) beam and pure 
gas targets internal to the HERA storage ring. The beam energy is 27.5 GeV 
and typical currents range from 40 mA at the start of a machine fill to 10-12 mA 
when the beam is dumped._The positrons are polarised after acceleration using 
the Sokolov-Ternov effect EH. This takes advantage of a small asymmetry in 
the spin-flip amplitude in the emission of synchrotron radiation. Polarisation 
builds up over a period of roughly half an hour to an equilibrium level between 
55% and 60%. The polarisation which is transverse to the beam momentum 
is rotated into the longitudinal direction in the HERMES straight section by 
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Table 4: Running conditions for the SLAC experiments. 



Experiment 


E s (GeV) 


PB 


Pr 




Target 


E142 


19-26 


0.36 


0.33 


0.33 


3 He 


E143 


9.7-29.1 


0.85 


0.7 (H) 


0.17 (H) 


NH 3 








0.25(D) 


0.25 (D) 


ND 3 


E154 


48.3 


0.82 


0.38 


0.55 


3 He 


E155 


48.3 


0.8 


0.83(H) 


0.23 


NH 3 








0.23(D) 


0.23 


LiD 



a series of magnets. It is rotated back to transverse downstream of the exper- 
iment. The beam polarisatiotLJs measured using Compton backscattering of 
circularly polarised laser lighted HERA is an e-p collider and while the proton 
beam is not used by HERMES, it nonetheless goes through the experiment, 
displaced horizontally by 72 cm from the positron beam. 

One of the strengths of HERMES is the purity of the targets. The dilution 
factor is 1 (the maximum of 1/3 for 3 He) so that the experimentally observed 
asymmetry is not quenched. Polarised gases (H, D, 3 He) flow into a windowless 
storage cell in the beam pipe. The cell has an elliptical cross-section (9.8 mm 
high x 29 mm wide) and is 400 mm long. It increases the areal density of 
the target by two orders of magnitude compared to a free atomic beam. Gas 
escapes out of the ends of the cell and is pumped away by a high speed dif- 
ferential pumping system. The 3 He target is polarised by optical pumpingEl. 
However, in contrast to the SLAC target, the 3 He is pumped directly so there 
are no extraneous atomic species to dilute the target. The 3 He nuclei are po- 
larised by spin exchange collisions with 3 S metastable 3 He atoms in a glass 
cell. The 3 S atoms are produced by a weak RF discharge and polarised by 
optical pumping with 1083 nm laser light. The target thickness is limited to 
10 15 nucleons/cm 2 by the requirement that the lifetime of the positron beam 
not be affected significantly by the target. The polarisation is measured by 
optical means and averaged 50% in 1995. The H and D targets are produced 
by an atomic beam source (ABS) based on the Stern-Gerlach effect to select 
specific spin states. The thickness of the H target was «7.10 13 atoms/cm 2 . 
The polarisation is measured by sampling the gas in the cell with a Breit-Rabi 
polarimeter which determines the fraction of atoms in each substate. Polari- 
sations as high as 92% were measured. However, recombination of atoms into 
molecules and depolarisation by collisions with the cell walls reduces this to 
«88% when averaged over all atoms in the cell. 

The HERMES detectorC-1 is a forward spectrometer with a large aperture 
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dipole magnet which provides a field of 1.3 T.m (see Fig. A horizontal iron 
plate shields the beams from the magnetic field. The spectrometer is therefore 
divided into two identical halves. Micro-strip gas chambers just downstream 
of the target exit window and drift chambers before and after the magnet mea- 
sure the scattering angle and momentum of charged particles. Proportional 
chambers inside the magnet track low momentum particles that do not reach 
the rear of the spectrometer. This is particularly valuable for pions from A 
decays. Particle identification is provided by a lead-glass calorimeter preceded 
by a pre-shower detector, a six-module transition radiation detector, and a 
threshold Cerenkov detector which is used primarily for pion identification. 
The Cerenkov detector has been upgraded to a dual radiator (C4Fio/aerogel) 
RICH during the 1997-98 shutdown. This will allow 7r-K-p separation over 
almost all of the kinematic range of the experiment (2-20 GeV). The large ac- 
ceptance of the spectrometer makes it possible to detect hadrons in coincidence 
with the scattered positron allowing HERMES to make the semi-inclusive mea- 
surements described in section |9| which are the cornerstone of the HERMES 
physics program. 

The running conditions at HERMES are summarized in Table 5. 
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Table 5: Running conditions for HERMES. 



Year 


E B (GeV) 


PB 


Pt 


Id 


Targi 


'95 


27.5 


0.5 


0.50 


1/3 


a Hc 


'96 


27.5 


0.55 


0.80 


1 


H 


'97 


27.5 


0.55 


0.88 


1 


H 


'98* 


27.5 


0.55* 


0.85* 


1 


D 


'99* 


27.5 


0.55* 


0.85* 


1 


D 



* Planned 



9 Semi-inclusive Deep Inelastic Scattering 

9.1 The Next Step towards Understanding the Spin of the Nucleon 

The experimental activity of the last decade has produced precise data for 
inclusive polarised DIS and can be summarized by the following statements: 

• The EMC result for the integral of g\ (x) has been confirmed. 

• The Bjorken sum rule has been verified to about 10%. 

• The integrals of g\{x), <?i(x), and gf(x) disagree with the Ellis- Jaffe 
predictions by up to several standard deviations. 

• The quark spins could contribute as little as 20% of the total nucleon 
spin. 

• The nucleon spin puzzle is not yet solved. 

The spin of the nucleon can be broken down into four components: 

S z = \ (AV + AS) + Ag + L z (44) 

where AV is the contribution from the valence quarks, AS comes from the sea 
quarks, Ag is the gluon polarisation, and L z is a possible contribution from the 
orbital angular momentum of the partons. Inclusive polarised DIS is sensitive 
mainly to (AV + AS). Since the quality of the inclusive data sets is now very 
good, further progress towards a solution to the spin problem must come from 
qualitatively different data, in particular hadrons produced in coincidence with 
the scattered lepton. Some information on Ag can be obtained from the Q 2 
evolution of g±(x,Q 2 ) but we have seen that these effects are small and that 
the current data are not precise enough to lead to a reliable result. 
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The strength of the semi- inclusive measurements is that the tyna.of hadron 
detected gives information on the flavour of the struck quark c3 ! E3. This is 
called flavour tagging. For example, a high energy n + coming from ae + -p 
collision was likely produced by the fragmentation of a u quark. This allows 
the spin distributions of particular quark flavours to be determined and terms 
in Eqn. |3| to be isolated. Some examples of this are discussed in later sections. 
Briefly: 

AV^: The valence quark distributions can be derived from data on the charge 



dependent polarisation asymmetry in pion production (Sec. 9.4) 



AS: The light sea quark distributions can be determined from the polarisation 
asymmetry of 7r~ production. 7r~ is a du bound state. Since the charge 
of the u is twice that of the d, it dominates the cross-section (remember 
the term) and this process is sensitive to the light sea. Information 
on the stra nge sea can be obtained from K~ production since it is a us 
state (Sec. |9.4| ). 

Ag: The photon-gluon fusion process dominates the production of charm par- 
ticles in DIS. Information on Ag can therefore be obtained from the 



polarisation asymmetry in J/tp and D° production (Sec. 9.6) 



L z : It has been proposed that the orbital angular momentum contribution to 
the nucleon's spin might be deduced from an azimuthal asymmetry in 
hadron production with a transversely polarised target. However, this is 
still controversial. 

9.2 Kinematics; the Current and Target Regions 

Two new variables are used to characterize hadrons as coming either from the 
struck quark [current region) or from the target remnant [target region): 

z = - 5 (45) 

z is the fraction of the energy transferred in the DIS process (v) carried by 
the hadron, and x p (x-Feynman) is a measure of the longitudinal momentum 
of the hadron in the 7* — p reference frame. Current quarks tend to produce 
hadrons at high z and positive Xf, while the target region is characterized by 
low z, negative- x p hadrons. We will see later that we are interested primarily 
in hadrons produced by current quarks. 

The question of how high in z is enough to ensure that the hadron comes 
from the current quark is not trivial. Surely z > 0.5 is safe but most hadrons 
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Figure 15: z-distribution of hadrons at HERMES in 1996. 



are produced below z = 0.5 (see Fig. [l5]) and lowering the cut on z is necessary 
to obtain a sufficient number of events. This problem was studied by Berger 
who defined criteria to ensure that a given hadron came from the current 
quark eZL He considered the rapidity: 

Vh = o lo M ~z? — ( 46 ) 



2 \E h -Ph,\\ 

and used the experimentally observed fact that the typical hadronic correlation 
length in a jet is 2 units of rapidity. In other words, it is highly unlikely that 
two jets separated by Ayh > 4 will contain hadrons coming from each other. 
Since the maximum range in yh which is kincmatically allowed for a jet is 

Y = log W 2 = log (V i-^) , (47) 

the condition Y > 4 translates to W > 7.4 GeV. Therefore, if we look along the 
direction of the virtual photon and apply Berger's criterion, we are reasonably 
assured that we are looking only at hadrons from current quarks. However, 
it is important to realise that this criterion need only be satisfied if we want 
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to consider all values of z. The condition on W can be relaxed if a cut on 
z is used. For example, for Y > 2 (i.e. W > 3 GeV) it is likely that the 
hadron which is carrying the most energy (called the leading hadron) comes 
from the current quark. Another useful guideline is that for z > 0.2, we need 
only require W > 4.8 GeV. 

In practice, the best way to be sure that contamination from the target 
region is not significant is to determine the quantity of interest with varying 
cuts on z and show that the result is stable. 



9.3 Factorisation and the Semi-inclusive Cross-Section 
The multiplicity of hadron h in the final state is given by: 

J_ d 3 cr h = X^e? gj{x,Q 2 ) D?(z,Q 2 ) 
a to t dxdQ 2 dz J2i e i Qi( x >Q 2 ) 



(48) 



where the sum is over quarks and anti-quarks separately, atot is the inclusive 
cross-section, is the charge of the struck quark, qi(x,Q 2 ) are the quark dis- 
tribution functions, and D^(z, Q 2 ) are the fragmentation functions. £>' l (z, Q 2 ) 
is the probability density that a quark i will produce hadron h with energy zv 
when it undergoes the process of fragmentation as it leaves the nucleonlia. In 
other words, Eqn. ^8] tells us that the mean number of hadrons h produced is a 
convolution of the probability to strike a quark qi and the likelihood that this 
quark will fragment into hadron h. In principle, the fragmentation functions 
depend on Q 2 as well as z but this is a relatively small effect and is ignored in 
most of what follows. Q 2 effects in fragmentation are studied as corrections to 
the final results. 

Implicit in Eqn. ^| is the concept of factorisation. This means that the 
fragmentation process is independent of the hard process that produced the 
outgoing quark, which in DIS is the absorption of the virtual photon by the 
struck quark. In practice, this means that in Eqn. ^| the fragmentation func- 
tions depend only on z. This leads to the remarkable fact that the fragmen- 
tation functions measured in e + — e~ collisions should h£.the same as those 
measured in DIS. This has been verified experimentally e3. Also assumed in 



Eqn. |48| is that the hadrons observed are from current quarks. If the latter 
condition is not met then factorisation might appear to be broken. 

Ultimately, the best way to be sure that factorisation holds is to plot an 
x-dependent quantity as a function of z and show that it is independent of this 
variable. The same can be done for z-dependent functions in terms of their 
x-dependence. 
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Figure 16: Favoured (D\) and unfavoured (D2) fragmentation functions from HERMES 50 
(filled circles) and EMC 51 (open circles). The EMC data have been evolved to the same Q 2 

as HERMES (2.35 GeV 2 ). 



Some simplifying assumptions can be used to reduce the number of terms 
in Eqn. For example, consider ir^ production on the proton. Even if 
only three quark flavours (u, d, and s) are included, we have 12 fragmentation 
functions. Using isospin symmetry and charge conjugation invariance, we can 
make the following simplifications: 

Di = Df = Df = Df = Df (Favoured FF) (49) 

D 2 = Df = Df = Dl = Df (Unfavoured FF) (50) 

£> 3 = Df = Df = Df = Df (Strange FF), (51) 

reducing the number of fragmentation functions to 3. Eqn. ^ states that the 
probability of a u or a d quark to produce a n + is the same as for a d or a u quark 
to produce a 7r - . Since these are valence quarks in the final state hadron, this 
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fragmentation function {D\) is called favoured. The processes referred to in 
Eqn. [5(] are more complicated and therefore D2 should be smaller than D±. 
D2 is called the unjaunured fragmentation function. Plots of D\(z) and D%(z) 
from the HERMES B and EMC El experiments are shown in Fig. [l6| where 
it can be seen that D\ is several times as large as D2 in the mid-z range. 
Finally D3 is the strange fragmentation function and should be similar in size 
to D 2 . Fragmentation (or hadronisation) is governed by long-range dynamics. 
It is therefore non-perturbative and we must use phenomenological models tq 
describe it. One such simple model is the independent fragmentation modeled 
which leads to: 

D 1 (z)~l + z- (& ' 

Eqn. [52] predicts D1—D2 at z — where the target region dominates and the 
hadronisation process is complicated enough that 7r _ 's should be produced 
in the same amount as 7r + 's. Eqn. |52| also predicts D 2 =0 at z = 1. This 
is intuitively correct since a hadron carrying all of the energy transferred in 
the reaction must contain the struck quark and the unfavoured fragmentation 
function should therefore vanish. Thia^behaviour is_seen in Fig. [It] which is a 
plot of D 2 {z)/D 1 (z) from HERMES E3 and EMCB. The data indicate that 
D 2 (z) / D\(z) does not drop off with z as quickly as expected from Eqn. p2[ , 
especially for z > 0.7. However, the error bars are large and the high z 
behaviour remains to be confirmed. 



9.4 Polarised Valence and Sea Quark Distributions 
Valence Quarks: 

The contribution to the nucleon spin from the valence quarks can be deter- 
mined, from the polarisation dependent yield of charged pions from two tar- 
gets e3. We will work through this example in some detail to show how an 
apparently complicated quantity can be written simply in terms of four quark 
distribution functions. 

First, write out the numerator of Eqn. ^8] for n + production on a H target 
when the spins of the beam and the target are opposite: 

Nft = ^u+DZ + +±d+Df +±u+Dt ^-d+Df+^s+Df +\s+Df . (53) 

Using Eqns. [l9[ and [51], we can re- write this as: 

= ^u+D l + ^d+D 2 + \u + D 2 + -d+D! + i.s+Z? 3 + ^ + D 3 . (54) 
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Figure 17: Plot of the ratio of the favoured and unfavoured fragmentation functions which 
qualitatively shows the behaviour predicted by Eqn. 

A similar equation holds for the production of 7r~ on H with the same relative 
spin orientations: 

Nf[ = ^u+D 2 + \d+D 1 + ^u+D 1 + \d + D 2 + ^s+D 3 + is+D 3 . (55) 



Taking the difference of Eqns. B4l and 55, the strange quark contributions drop 



out and the term {D\ — D 2 ) factors out: 

- JVff = \{u + Q+)(i?i - AO - \{d+ d+)(D 1 D 2 ). (56) 



Since {q + — q + ) is just the valence distribution for quark flavour q, Eqn. 56 
can finally be written as: 

Nfi - Nf~ = l(4u+ - d+)(Ua - D 2 ). (57) 

The sea quark contribution has been cancelled out by taking the difference 
of 7r + and 7r~ production and we are left with an expression for the valence 
quarks alone. A similar expression holds for the other spin direction, which is 
sensitive to u~ and d~: 

Af T + - ATff = I(4u" - - D a ). (58) 
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Taking the difference (sum) of Eqns. ^7] and [38], we get an expression which 
depends on the polarised (unpolarised) valence quark distribution functions: 

(7Vf x + - Nf[ ) - (7Vf T + - Nf ) = ^(4A Wl) - Ad v )(D! - D 2 ), (59) 

(iVf/ - Nf[ ) + (iVf T + - Nf ) = ^(Au v - d v )(D 1 - D a ). (60) 

Finally, the fragmentation functions cancel in the ratio of Eqns. [59| and |60| and 
we get: 

pjzisnz^zm^ = (61) 

(iVf^ - ) + (A^+ - N$ ) L 4u„ - d v 

The unpolarised quark distribution functions in the denominator have been 
determined from unpolarised DIS experiments so that we are left with only 
two unknowns in Eqn. [u]. A similar expression can be written for a deuterium 
target: 

pLZ K )-(<-*?,-) [deuteto „, _ Au, + Ad 

and we can now solve for Au v and Ad v separately. This is a powerful tool 
to isolate the contribution to the nucleon spin from the valence quarks. Note 
that we must be able to identify pions for this technique to work. However, 
the same formalism can be applied to all hadrons at the expense of extra 
model dependence if Monte Carlo techniques are used to calculate corrections 
for the fact that the type of final state hadron is not determined. Note also 
that data on two targets are needed to determine the polarised valence quark 
distributions using this technique. 

Sea Quarks: 

Different combinations of identified hadron yields can be used to isolate dif- 
ferent quark distributions. In particular since tt~ is a ud state and that the 
coupling in DIS scattering is four times larger for the u quark than for the d 
quark, polarisation dependent 7r~ production should be sensitive to the po- 
larised u distribution. Furthermore, since K~ is a us state, i±_is all sea and 
should have increased sensitivity, especially to the strange seae£l. Expressions 



for 7r and K production based on Eqn. 48 can be derived in a similar way 
to the previous section. 
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HERMES Semi-Inclusive Hadron Asymmetry on the Proton 
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Figure 18: Single hadron asymmetries from SMC 52 and HERMES 53 for a H target. 

Once high-statistics data sets exist on the production of -k^ and K ± on 
more than one target, a global fitting procedure can be used to determine all 
polarised quark distributions at once (see next section). It is hoped that such 
data sets will be available in the next two years. 

9.5 Single Hadron Asymmetries and Polarised Quark Distributions 
Single Hadron Asymmetries 

Semi-inclusive data are usually presented in the form of single hadron asym- 
metries defined as follows: 

A h = — U (63) 

N u + N n 
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HERMES Semi-Inclusive Pion Asymmetry on the Proton 
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Figure 19: Single pion asymmetries from HERMES 53 for a H target. The data are compared 
to quark distributions from a leading order fit by Gehrmann and Stirling 30 . 



Data on semi-inclusive processes exist only from HERMES and SMC, with 
identified pions only available at HERMES. The acceptance of the SLAC spec- 
trometers makes detection of a hadron in coincidence essentially impossible. 
SMC has data for both H and D, while HERMES currently only has data on 
H and 3 He (D will be run in 1998-99). The charged hadron single spin-asym- 
metries for the proton are shown in Fig. |l| for SMCE3 and HERMES EE The 
error bars are comparable. However only one third of the HERMES data are 
included in this plot. As in the inclusive case, the SMC data go to lower x 
because of the higher beam energy. The data quality is already quite good and 
will be improved when HERMES finishes the analysis of the 1997 data. The 
asymmetries are positive and rather large, especially at x > 0.1. 

Pion asymmetries from HERMES @ are shown in Fig. Il9| where they are 
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Figure 20: Single hadron asymmetries from HERMES 54 for a 3 He target. 



compared to the fit of quark distribution functions by Gehrmann and Stir- 
ling l3. Note that here too, only one third of the data have been analysed. 

Finally, Fig. ^0| shows the hadron asymmetries measured by HERMES l-S 
on a 3 He target in 1995. Since this was the commissioning year for the experi- 
ment, the statistics are quite low and the error bars are correspondingly large. 
However, it is clear that the asymmetries arc much smaller for the neutron 
than for the proton as in the inclusive case. The data set for the neutron will 
be improved significantly in 1998-99 when HERMES runs with a deuterium 
target. 



Polarised Quark Distributions; the Purity Method 

The best way to determine the polarised quark distribution functions has been 
to perform a global fit to the inclusive data. The same can be done includ- 
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ing the semi-inclusive data when a sufficiently large data base has been ac- 
cumulated. We can write the single hadron asymmetries in terms of quark 
distribution functions and fragmentation functions: 

Ah I N \\ N \\ EA ^Q 2 ) Dl^Q 2 ) (1 , R(x Q2)) m 

Al -DN? l+ N^- Y.A «{x,Q>) m*,Q 2 ) (1 + i?(x ' 0)) ' (64) 

where the sum is over quarks and anti-quarks separately and R = (Jl/^t as 
seen earlier. The (1 + R) term is needed to be consistent with the extraction of 
the unpolarised quark distribution functions using F 2 which does not assume 
R = 0. We can define functions P/ 1 , called Purities, as the probability that 
a hadron h was produced by the fragmentation of a quark qi that was struck 
in the DIS process E3. These are similar to the fragmentation functions but 
should not be confused with them. The purities are written as follows: 

P h (x Q 2 z) = e > 2 g^,Q 2 ) jL D H*)dz 
We can now re- write Eqn. in terms of purities and quark polarisations: 

a1( x ) = J2p«(x) ^rv C 1 + (66) 

q q[ ~ X) 

The purities can be calculated by Monte Carlo and hence can also contain 
information about the acceptance of the spectrometer and detector effects. 
Purities can also be defined for inclusive asymmetries as follows: 

„ e 2 q(x,Q 2 ) 

and give the probability that the inclusive scattering came from quark q. If 
one combines inclusive and semi-inclusive measurements, Eqn. |6^ becomes a 
matrix equation and we can determine the polarised quark distribution func- 
tions by inverting the purity matrix. For example, consider the asymmetries 
for inclusive positrons (Af ), and charged pion production (A\ and A\ ): 



/ 


At N 






Af 


H 


V 








(68) 



The (1+R) term is now included in the matrix P. More data can be included to 
get a larger matrix equation and determine more quark distribution functions. 
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Figure 21: Polarised quark distributions from a fit to SMC data 52 and using the purity 
method for HERMES 1996 data 55 . Distributions are plotted for valence u and d quarks 
as well as for the sea q under the assumption that the polarisation of the u sea , d sea , and s 

quarks is the same. 
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Care must be taken when inverting the purity matrix because the data are 
correlated (e.g. the semi-inclusive data are a sub-set of the inclusive data) 
and the uncertainties must be calculated correctly. For more details, see the 
contribution to this workshop by M.-A. Funk. 

Results from SMC@ and HERMES (1996 only) HI are shown in Fig. ^ for 
the polarised valence quark distributions as well as for the sea distribution q. 
It is not possible with the current data to separate the contribution from each 
type of sea quark. New data to be taken in the next two years at HERMES 
with the RICH detector and with two targets (D and H) will greatly improve 
the situation. 

9.6 Charm Production and the Gluon Contribution 

While data collected to date on inclusive and semi-inclusive polarised DIS 
are sensitive to the spins of the quarks, we saw earlier that it is now crucial 
to determine if the gluons are polarised. In order to be sensitive to gluon 
polarisation, the normal DIS diagrams must be suppressed. One way to do 
this is by considering charm production. This has proven to be useful in the 
unpolarised case E3. There should not be a large charm component in the 
nucleon wavefunction (intrinsic charm) so it is not very probable to strike a 
charm quark directly. Furthermore, charm production in the fragmentation 
of light quarks is suppressed by the mass of the charm quark. As a result, 
charm production in DIS is dominated by photon-gluon fusion which is shown 
diagrammatically in Fig. |2^. In this process, the virtual photon and the gluon 
interact through the intermediary of a cc pair. The charm quarks can manifest 
themselves as charmed mesons in the final state (open charm production: D- 
mesons), or as 3/tp mesons (hidden charm). In the latter case, a second gluon 
must be emitted by the cc system in order to get a colour singlet final state. 
If this gluon is soft, the interpretation of the process becomes more difficult. 
Therefore, conditions must be placed on the kinematics of the reaction to 
ensure that we have inelastic J/^'s in the final state. The asymmetry in charm 
production for the two relative spin states of the beam and the target can be 
written as: 

Ag 

All{z,pt) = &LL— (69) 

where px is the transverse momentum of the charm meson, and (111 is the 
asymmetry for the elementary hard scattering process (i.e. 7 — c). The asym- 
metry depends on the production mechanism. It has-been calculated for J ftp 
production in the colour singlet modelc3 by Guilletcl and is « 0.5 for HER- 
MES kinematics, qll has also been estimated in the colour octet model which 



41 



e 




Figure 22: Diagram of the photon-gluon fusion process. The photon and the gluon couple 
through the intermediary of a cc-pair. 
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Figure 23: J/i/) mass peak in the e + -e and [i^-fJ, decay channels at HERMES. The 
efficiency for /i identification will be significantly improved for f998. 
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might be more applicable at high z. In this case oll=1 so that the asymmetry 
is likely to be large in either model. J/^'s have been clearly identified at HER- 
MES (Fig. ||) and the projected_accuracy on Ag/g is 0.52 after the '98-'99 
running period and 0.38 by 2001 Ea. This will be reduced if improvements to 
the target can be made. While this accuracy is limited, HERMES should be 
able to make a qualitative statement on whether Ag is non-zero. 

An advantage of open-charm production is the larger cross-section com- 
pared to J/ip production. This is particularly important at HERMES where 
the cm. energy is barely above threshold. On the other hand, the identifica- 
tion of D-mesons is more difficult than J/?/>, especially without K/n separation 
in the hadron PID. 



10 The Future of the Field 

HERMES is the only running polarised DIS experiment for the moment. SLAC 
has approved an extension to E155 to collect high statistics data with a trans- 
versely polarised target to improve the data on g 2 . This experiment should run 
in 1999. An experiment called COMPASS is being prepared at CERN to mea- 
sure semi- inclusive DIS, in particular open charm production at higher beam 
energy than HERMES to study the gluon polarisation. The capability to have 
polarised protons in RHIC is being included in the accelerator design. The 
physics to be addressed in high energy p-p collisions is described briefly below 
and in another contribution to these proceedings Ea. Finally, DESY is studying 
the possibility to accelerate polarised protons in HERA to study polarised DIS 
at very small values of x. 

10.1 HERMES Upgrades 

HERMES has embarked on an upgrade program for the 1998 run and beyond 
to include kaon identification and improve the efficiency for the detection of 
charm mesons e3e30. A schematic of the HERMES upgrades is shown in 
Fig. §|. 

Hadron identification will be greatly improved by replacingJdie threshold 
Cerenkov detector by a Ring Imaging Cerenkov detector (RICH)EJ. The RICH 
will use a dual radiator system: C4F 10 gas with an index of refraction of 1.0014 
and clear aerogel with n=1.03. The Cerenkov photons will be detected by an 
array of 3/4" phototubes (roughly 2000 per half). The RICH will provide 
identification of pions, kaons, and protons (and p) over most of the kinematic 
range of the experiment (2-20 GeV). This will allow HERMES to measure all 
hadron asymmetries and thus determine the contribution of each quark species 
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Figure 24: Schematic representation of the upgrades to the HERMES spectrometer. The 
RICH, the iron wall for muon identification, and the FQS to measure the kinematics of 
electrons from low Q 2 events will all be ready for the 1998 running period. The limited 
streamer tubes to increase the acceptance for fi + -fi~ pairs and most of the silicon-strip 
detectors in the target chamber will not be ready until 1999. 



44 



to the nucleoli's spin. Furthermore, kaon identification will allow the detection 
of D-mesons through the K-ir decay channel. 

Further JLmprovements to the charm detection efficiency will be made on 
three frontsLJ. First, an iron wall to absorb hadrons has been erected behind 
the calorimeter followed by a scintillator hodoscope. This system will improve 
the muon identification and double the number of J /tp particles detected by in- 
cluding the /i + /i _ decay channel. Better muon identification will also help for 
semi-leptonic D-meson decays. The second upgrade related to charm is a For- 
ward Quadrupole Spectrometer (FQS) which will be installed in the first two 
machine quadrupoles downstream of the experiment. This device will detect 
electrons from low-Q 2 events, which are scattered down the beampipe. This 
will allow full kinematic reconstruction of these events and help separate elastic 
(or diffractive) J /tp production from events which can be interpreted in terms 
of Ag. Finally, the acceptance for J/ip — > will be increased by extending 

the tracking in front of the magnet and instrumenting the corresponding region 
behind the magnet steel which will act as a muon filter. The enlarged /i + ji~ 
acceptance will greatly increase the 3/tjj sample since most of these particles 
are produced with very low energy and hence with a large opening angle for 
the decay muons. 

HERMES is also preparing tap arrays of silicon-strip detectors to be lo- 
cated inside the target chamber EH: two wheel-shaped detectors between the 
storage cell and the exit window (A- wheels), and recoil detectors parallel to the 
storage cell. The recoil detector has many functions. It can detect slow recoil 
protons in e-d scattering to tag e-n events and eliminate nuclear corrections. 
The detection of a recoil proton in J/-0-production can be used to eliminate 
elastic events and also to help reconstruct decays where only one of the leptons 
is in the spectrometer acceptance. As their name implies, the A wheels will 
be used to study A production and the subsequent decay to p + ir~ (A) or 
p + 7r + (A). The idea is to take advantage of the fact that the A-decay is 
self- analysing making it possible to measure the polarisation of the A's. This 
polarisation can be written as: 

P A = (P B D + P T — ) (70) 

It IT 1 

where Pb and Pt are the beam and target polarisations respectively, D is a 
depolarisation factor and Au/u and Au A /u A are the polarisations of the it- 
quark in the proton and the A. With an unpolarised target, Au A /u A can be 
determined and the spin transfer from the it-quark in the proton to the A can 
be measured. If we assume that this quantity does not depend on the target 
orientation, we can then determine the transverse polarisation of it-quarks 
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in the proton by measuring Pa with a transversely polarised target. While 
the current acceptance of the HERMES spectrometer is not well matched to 
A detection, the A-wheels will substantially increase the acceptance and also 
extend the range for negative values of Xf- 



10.2 COMPASS at CERN 

The COMPASS experiment at CERN is being designed to study_semi-inclusive 
processes at a higher energy than HERMES (100 GeV muons)l£3. The goal is 
to measure Ag/g to an accuracy of 0.11 after 1.5 years of data taking. The 
main advantage over HERMES is the much higher charm production cross- 
section due to the higher beam energy. The principal channels to be studied 
are: 

A two stage spectrometer with very large acceptance is being built. Good 
hadron identification will be provided by two RICH detectors. The targets will 
be similar to those used at SMC and E155: NH3 (85% polarisation) and LiD 
(50% polarisation). COMPASS will also make the semi- inclusive measurements 
described in Sec. 0. The experiment is expected to start data taking in 2000. 



10.3 Polarised beams at RHIC 

The capability to study collisions of 250 GeV polarised protons at RHIC is 
being prepared Ej. A more complete description _af this program is given in 
the contribution to these proceedings by J. HarrisEj and the reader is referred 
to this paper. Briefly, the polarisation of the sea quarks can be measured 
using the Drell-Yan process (qq — > while sensitivity to the gluon po- 

larisation can be obtained by studying the production of isolated photons (the 
so-called direct photon process: qg — > 73), jet production, or i/ip production. 
The direct photon experiment is challenging because of the large background 
of 7's from 7T° decays. While it is possible to isolate gg-processes by using 
the characteristics of jets (e.g. the invariant mass of the 2-jet system), the 
interpretation is more difficult than for the direct photon process. RHIC is 
scheduled to produce heavy ion collisions in 1999 and it is expected that some 
time could be set aside for polarised protons as early as 2000. 



10.4 Polarised Protons at HERA 



In the long term («2005), . 
polarised protons in HERA I 
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pY is studying the possibility of accelerating 
1 which would allow the extension of the data 



set on gi down to very low x. One of the problems associated with this program 
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is the small dcpolarisation factor D due to the kinematics of the scattering. 
This means that the polarisation of the virtual photons is small and very high 
luminosity is required to compensate. Furthermore, getting polarised protons 
through the accelerator chain at DESY is not trivial and will involve the use of 
partial Siberian snakes in PETRA and as many as four full snakes in HERA. 
Nevertheless, information on g\ at low x is crucial to our understanding of 
how to do the extrapolation to x = 0. HERA would extend the minimum x 
from about 5.1CP 3 down to a few times 10 -5 . Furthermore, asymmetries in 
the production of jets are expected to be large and will give information on 
Ag. 

10.5 Real photon experiments 

Finally, experiments on charm production using real photons have been pro- 
posed both at SLAC (E156) and at DESY (APOLLON). However, neither of 
these has been approved and they are on hold until further notice. 

11 Summary 

A pedagogical overview of the field of nucleon spin physics has been presented 
in these lectures. Ongoing technological advances in polarised target and beam 
production have made possible a series of experiments over the last 10 years 
to study how the spin of the nucleon is related to the spins of its constituents. 
There now exists a large body of data on inclusive polarised deep inelastic 
scattering from CERN, SLAC, and DESY. These data indicate that only a 
small fraction of the nucleon's spin comes from the spins of the quarks. While 
the Bjorken sum rule has been verified to about 10% accuracy, the data on the 
integrals of gi (x) for the proton and the neutron disagree with the predictions 
of the Ellis- Jaffe sum rules by one to several sigma. Further progress in this 
field will be based on semi-inclusive measurements where the flavour of the 
struck quark can be tagged by the type of hadron produced in the final state. 
This will allow the contributions from valence and sea quarks to be isolated. 
The production of charmed mesons will give information on the polarisation of 
the gluons. The HERMES experiment which will run with a deuterium target 
in 1998-99 will produce complete data sets on semi-inclusive pion and kaon 
production for the proton and the neutron. HERMES will also study charm 
production and should be able to make a qualitative statement on whether 
the gluons are polarised. A precise determination of Ag must await the new 
experiments being planned at CERN (COMPASS) and at Brookhavcn (Spin- 
RHIC). The next few years promise to be very exciting in this field and we 
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might finally find a solution to the "nucleon spin puzzle" . 
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